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The Bolztmann echo (BE) is a measure of irreversibility and sensitivity to pertur- 
bations for non-isolated systems. Recently, different regimes of this quantity were 
described for chaotic systems. There is a perturbative regime where the BE decays 
with a rate given by the sum of a term depending on the accuracy with which the 
system is time-reversed and a term depending on the coupling between the system 
and the environment. In addition, a parameter independent regime, characterised by 
the classical Lyapunov exponent, is expected. In this paper we study the behaviour 
of the BE in hyperbolic maps that are in contact with different environments. We 
analyse the emergence of the different regimes and show that the behaviour of the 
decay rate of the BE is strongly dependent on the type of environment. 
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1. Introduction 

In quantum mechanics there is no "exponential separation" of initial conditions 
due to chaotic motion because evolution is - in principle- unitary. Peres (1984) 
proposed, as an alternative, to study the stability of quantum motion due to per- 
turbations in the Hamiltionian. As a consequence, the Loschmidt Eco (LE) (Peres 
1984; Jalabert & Pastawski 2001; Jacquod et al. 2001; see two reviews: Gorin et al. 
2006 and Petitjean & Jacquod 2009) 



M{t) = (^o|e'^^*/''e-*^*/''|V'o) 



(1.1) 



was introduced with the purpose of characterising the sensitivity and irreversibility 
arising from the chaotic nature of quantum systems. The parameter S denotes 
perturbation strength. Equation (jl.ip has a dual interpretation. On the one hand, 
it can be interpreted as how close a state remains to itself evolving under slightly 
different Hamiltonians. On the other hand, it measures the sensitivity of a system 
to imperfect time inversion, i.e. evolve forward in time under Fl and then invert 
time and evolve backward with iJs (supposing that the time inversion operation is 
not perfect). 
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Depending on the nature of the underlying dynamics, the LE can exhibit quahta- 
tively different behaviour and it thus can be used to characterise quantum chaotic 
systems. Moreover a number of time-reversal experiments have been performed 
(Hahn 1950; Rhim et at 1970; Zhang et al. 1992; Pastawski et al. 2000), and therein 
lies the importance of the LE. In addition, the LE (which in quantum information 
is known as fidelity) can be efficiently measured in quantum information systems, 
i.e. its measurement scales only polynomially with the system size (Emerson et al. 
2002). 

An important fact to remark is that quantum systems cannot be isolated easily. 
Most of the times, there is an environment acting upon the system. This interaction 
is most likely unknown and its effects may be uncontrollable. The Boltzmann echo 
(BE) was introduced (Petitjean & Jacquod 2006) as a generalisation of the LE 
to take into account the fact that quantum systems are not isolated. The idea is 
to consider the evolution of a system s with a Hamiltonian Hg which is coupled 
to a an environment e whose evolution is given by Hf,- We suppose the evolution 
of the environment e is unknown and are therefore uncontrollable, so we trace 
out the environment degrees of freedom. Given a separable initial state, such as 
Pa = p'o'' P^q \ where we take Po*' — |'0o)('0o|, the BE is defined as the partial 
fidelity 

MBit) = {{i:o\^v,[e-'"'''^e-'"f'l^poe'"^'l^e'"''^']\^o)), (1-2) 
where i/f, and Hf are given by 

Hf = H.,®l^+ls®H^ + Uf (1.3) 

Hb = -{Hs + S,) ® le + I, ® ~{H, + Se) + C/b, (1.4) 

and represent the forward and backward Hamiltonian respectively. Equation (jl.2p 
can be explained as follows. First take an initial state po and evolve it forward up 
to time time t with Hamiltonian Hf. Then, invert time evolution and evolve with 
Hamiltonian i?f,. The imperfection in the inverting process is represented by: Eg 
for the system. Eg for the environment. The terms Uf,Ub represent forward and 
backward interaction between system and environment (for simplicity throughout 
this work we consider Uf = ~Ub). Finally, the evolution of the system and the BE is 
obtained by performing a partial trace over the environment degrees of freedom and 
computing the overlap. Tracing out the environment makes the effective evolution 
of the system non- unitary producing decoherence (Zurek 2003). An average over 
initial states of the environment pg'^'' is necessary (represented with big brackets in 
equation (jl.2p ) because we have no control over its degrees of freedom. 

In the work of Petitjean & Jacquod (2006) the BE was studied semiclassically 
for two interacting - classically chaotic - sub-systems. One of them was used as 
system and the other as an environment. They found three different regimes for 
the BE as function of time: parabolic or Gaussian for very short times; exponential 
for intermediate, followed by a saturation depending on the effective Hilbert space 
size. Here we focus on the exponential regime and specifically on the dependence 
of the decay rate on the perturbation and environment parameters. The authors 
show (Petitjean & Jacquod 2006, see also Petitjean & Jacquod 2009) that in the 
Fermi golden rule (FGR) regime (small perturbation and weak coupling with the 
environment) the decay rate of the BE results from the sum of the decay rates 



Article submitted to Royal Society 



Irreversibility in quantum maps with decoherence 



3 



of the LE due to imperfect time inversion (by definition the BE in the limit no 
decoherence is just the LE), and the contribution due to the interaction Uf, Ub 
with the environment, F — + F/ + F^. Henceforth, we cah this the sum law. 
Moreover, for chaotic systems they find that in the hmit of strong environment 
couphng or large perturbation the decay rate is perturbation independent and is 
given by the classical Lyapunov exponent. 

In the present contribution, we study the BE for quantum maps on the torus 
that are classically chaotic. Quantum maps are very simple models that have all 
the main features of chaotic systems and are ideal for numerical studies. Our goal is 
to understand the behaviour of the BE under the action of different environments. 
For this reason we have computed the decay of the BE for a wide range of the 
parameters that control the perturbation of the system and the interaction with 
the environment. We find that a sum law for the decay rate of the BE exists. It 
can be expressed as the sum of the decay rates of the LE and the purity of the 
system, but it is fulfilled only partially, depending on the decoherence model. The 
decoherence models that we present can be written as a convolution with a kernel. 
It is for the cases where the kernels have polynomially decaying tails -models with 
somewhat large correlations in phase space- when the sum law is best achieved. In 
addition, the oscillations of the decay rate of the LE, found in e.g Wang (2004), 
Andersen (2006) and Ares & Wisniacki (2009), are damped completely in the limit 
of strong decoherence. However, the decoherence (and perturbation) independent 
decay rate saturation at the classical Lyapunov exponent is not present for all 
decoherence models. 

The paper is organised as follows. In SjHHwe describe the quantum kicked maps 
on the torus, the systems used for our studies. Then in, we introduce our 

model of open maps using translations in phase space and the Kraus operator sum 
form. The main part of this contribution is which is devoted to the numerical 
calculations and presentation of the results. Finally, in 21 we summarise our work 
and results. 

2. The system 

(a) Quantum 'kicked' maps 

Classical maps generally arise from the discretisation of a differential equation 
of the motion - like e.g. a Poincare surface of section. Nevertheless one can build 
abstract maps that do not necessarily relate to a differential equation but that can 
however provide insight into the properties of chaotic dynamics - e.g the baker's map 
or the cat map. Like classical maps, quantum maps are usually simple operators 
with all typical properties of quantum chaotic systems like level spacing statistics. 
In addition, there exist efficient quantum algorithms for some quantum maps (e.g. 
Goergeot & Shepelyansky 2001; Levy et al. 2003). As the Hilbert space grows 
exponentially with the number of qubits, one could reach the semiclassical limit 
with a relatively small number of qubits. For this reason they are ideal testbeds for 
current quantum computers in one of their possible uses: quantum simulators (see 
Schack 2006). 

The systems we consider are quantum maps on the 2-torus. Periodic boundary 
conditions imply that Hilbert space has finite dimension N and the effective Planck 
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constant is h= \/2'kN . This means that the semiclassical hmit is reached as N 

=0 



00. Position and momentum bases are discrete sets {g^ = i/^lilo^ ^'^'^ 



i/A^}^Q^ related by the discrete Fourier transform 

1 



= ^e-(^'^'^)^''. (2.1) 



For practical purposes we will consider maps which can be expressed as two 
shears -linear or non-linear- 



dV[q) 

1' = 



^ ^ 4p') ("lodl). (2.2) 



These maps can be quantised and the associated unitary map can be written as a 
product of two 'kicks' 

^ ^ ^i2TrNTip)^-i2^NV{q) ^ (^2.3) 

These types of map usually arise from Hamiltonians with periodic delta-kicks, like 
the kicked rotator (Chirikov et al. 1988) or the kicked Harper Hamiltonian (Leboeuf 
et al. 1990) . One of the advantages of implementing these types of maps numerically 
is due to the possibility of using the fast Fourier transform. 



(6) Open quantum maps 

A system with an evolution given by a map U might interact with another sys- 
tem acting as environment. If the dynamics of the environment cannot be accessed 
or controlled then the usual procedure is to trace out the degrees of freedom of the 
environment. Tracing out the environment translates into a loss information about 
the evolution, hence the word open - we picture information flowing out of the 
system. It is this loss of information the cause of decoherence - and subsequent loss 
of quantumness (Zurek 2003). For a Markovian environment and in the weak cou- 
pling limit, this is given by a completely positive-trace preserving map of density 
matrices into density matrices - sometimes called superoperator -which generally 
can be written in Kraus operator sum form (Kraus 1983) 

= (2.4) 

i 

where trace preservation is assured by K^Ki — / (/ is the identity)!^ Therefore 
the action of the environment is coded into the Kraus operators Ki, in analogy 
with the Lindblad master equation (Lindblad 1979) where the action of the en- 
vironment is given by the Lindblad operators. Different Kraus operators will give 
different types of environments. Rather than modelling the environment through 
the Lindblad operators and solving the master equation, here we directly model the 
effect of the environment on the density matrix of the system by 

N-l 

pt = D,(pt_i) = ^ Ce{q,p)TgpPt-iT^p, (2.5) 

p,q=0 

t Throughout this contribution the 'time' t is a discrete time variable which implies the number 
of times a map (or a superoperator) has been applied. 
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where T^p are the translation operators on the torus, Cc{q,p) is a function of q and 
p and e quantifies the strength of the system-environment couphng. Even though 
position and momentum operators with canonical commutation rules are not de- 
fined on the torus, translations can be defined as cyclic shifts over the bases ele- 
ments (Schwinger 1960). Since T^p are unitary, trace preservation in equation equa- 
tion (|2.5p requires that J2q p Cei q,p) — 1- The action decoherence superoperator 
Df introduced by equation (|2.5p has a simple interpretation: it implements every 
possible translation in phase space with probability Ce{q,p). This effect is clear in 
the Wigner function representation. Let W{q,p) be the discrete Wigner function 
(see e.g. Bianucci et al. 2002) of a density matrix p then equation equation (12. Sp 
can be re- written as a convolution with Cg {q, p) 

W^t(g,P) = ^c,(<z,p)l¥t_i(g-g,p-P). (2.6) 

q,p 

This is an incoherent sum of slightly displaced Wigner functions. Any fast oscillating 
term present in the state represented by W{q,p) will be eventually washed out, 
depending on the form of Ce{q,p). 

For simplicity we suppose that the complete evolution of the quantum map and 
the decoherent part take place in two steps: first the unitary map U followed by 
the decoherence term of equation (|2.5p 

pt = D,{Upt^iU^). (2.7) 

This is an approximation that works exactly in some cases, e.g. a billiard that has 
elastic collisions on the walls and diffusion in the free evolution between collisions. 
This kind of two-step model has been used to study quantum to classical corre- 
spondence and the emergence of classical properties from the quantum dynamics 
(Nonnenmacher 2003; Garcia-Mata & Saraceno 2004) 

The effect of decoherence can be characterized by using the purity 

P{t) - tr(p2), (2.8) 

were pt is the reduced density matrix of the system. The purity measures the relative 
weight of the non-diagonal matrix elements. It is a basis dependent measure that 
can be used to quantify the amount entanglement between two parties. If P{t) = 1, 
it means that the global system can be factorized into two separate systems and 
there is no entanglement. On the contrary, if the purity of the reduced density 
matrix is minimum (completely mixed state), then the entanglement is maximal. 
In the case of an N dimensional system P{t) ~ 1/N for a completely mixed state 
(maximally entangled with the environment). As a function of time, after an initial 
short transient, the purity decays exponentially. For long times it saturates to a 
minimum value given by fi,/(27r). 



3. Numerical results 

For our numerical calculations we use the cat map perturbed in position and mo- 
mentum with a smooth non-linear shear 

P', ^ P + ^r o";'wo"'^.^ (-odl), (3.1) 
q' = q + bp' - 2Trksin{2'!rp') ^ ^' ^ ^ 
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with a, b integers. This map is uniformly hyperbohc and fuhy chaotic. Fo fc ^ 1 the 
largest Lyapunov exponent given by A « ln((2 + a6+ y/ab(4: + ab))/2)/2. According 
to equation (j2.3l) the quantum version of equation (I3.ip is 

jj^ _ g27ri(-P^/(2Af)-fccos(27rP/iV)) g27ri(Q^/(2Af) + fe cos(27rQ/Ar)) 2) 

where P, Q = 0, ...,N— 1. All the arithmetic peculiarities of the cat map, which 
account for the non-generic spectral statistics are destroyed for fc 7^ ( Basilio 
de Matos & Ozorio de Almeida 1995; Keating & Mezzadri 2000). We can rewrite 
equation (|1.2I) for the BE for our open map as the overlap between two states 
evolving forward in time - with slightly different maps plus decoherence - as 

MbW -Tr[p-tpt] (3.3) 

where 

Pt = D,{UkPt-iUl), (3.4) 

Pt = D,{Uk'Pt-iUl,), (3.5) 

where k, k' are the perturbation strength of the cat map. We measure the pertur- 
bation of one map with respect to the other by the parameter 

Y.= \k'-k\. (3.6) 

For a chaotic system, after an initial transient the BE decays exponentially (Petit- 
jean & Jacquod 2009). Here we focus on the decay rate F as a function of E and e 
for the exponential decay regime. In the limit e — >■ we have F — Fs, where Fs is 
the decay rate of the LE. In the limit S — > the BE as defined in equation p.3p is 
equal to the purity, so decay rate is given by the decay rate of the purity F^ . 

We explore the behaviour of F for three decoherence models and a wide range of 
values of e and S. We analyse the parameter domain of validity of the sum law (now 
F = Fs + Ff) for these models. The different models of decoherence we consider are 
implemented simply by changing the coefficients Cc{q,p) in equation (j2.5p . Like for 
the LE, to extract the decay rate F an average over an ensemble of initial states 
needs to be performed. For the averages we used tt-s = 10 randomly chosen coherent 
states. 



(a) Gaussian diffusion 

The first model we have considered was introduced in the work of Garcia-Mata 
et al. (2003) to model diffusion in a quantum map. We take a periodic sum of 
Gaussians - to fit the boundary conditions of the 2-torus - 

1 

j,k——x 

where x is large enough (typically of order 10-15) so that the tails of the furthermost 
Gaussians can be neglected and A is the normalisation factor {q,p = 0, . . . , N ~ 1). 
We call this model Gaussian diffusion model (GDM). The GDM can be interpreted 



{q - jNf + {p- kNf 



2( 



(3.7) 
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0.005 0.01 0.015 

e 

Figure 1. (a) Decay rate F of the LB as a function of the rescaled strength of the pertur- 
bation E/ft for a GDM environment. The map is the quantum version of the perturbed 
cat [equation with a — b = 2. Averages were done over — 10 initial states. Other 

parameters are: k = 0.001, iV = 800. (x) e = (LE), (□) e = 0.003, (0) e = 0.0035, (A) 
e = 0.004, (O) e = 0.005, (V) e = 0.01. The horizontal dashed (in (a) and (b)) lines cor- 
respond to the Lyapunov exponents of the corresponding map A = ln[3 -f- 2\/2] ~ 1.76275; 
(b) The decay decay rate Fg of the purity as a function of the perturbation parameter e. 
The points correspond to the initial values of the curves in (a), (inset) Decay rate F — Fe 
as a function of the rescaled strength of the perturbation E/ft. 

as a smoothing or coarse graining of the unitary evolution: with Gaussian weight 
the state is displaced all over a region of size of order e. As a consequence, the 
interference terms get washed out, while the remaining classical part is diffused. As 
stated before, in the continuous limit equation (j2.6D is a convolution of the Wigner 
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function with a kernel Ce{q,p). For the GDM it can be related to the solution of the 
heat equation with diffusion constant given by (e/27r)^ (Zurek & Paz 1994; Strunz 
& Percival 1998; Carvalho et al. 2004; Wisniacki & Toscano 2009). 

In figure [Ha) we show the decay rate F of the BE as function of perturbation 
parameter S for the perturbed cat map a — b — 2, N — 800 and k — 0.01 in the 
presence of GDM for distinct values of e. The Lyapunov exponent A = ln[3 + 2^/2] 
is marked by a dashed line. For e = (red x symbol) we recover the decay rate of 
the LE: for small E we get the characteristic quadratic behaviour for small pertur- 
bation - Fermi golden rule regime; for larger values of E we get a non-universal - 
perturbation dependent ~ oscillatory behaviour which has also been observed in the 
work of Wang et al. (2004), Ares & Wisniacki (2009) and Casabone et al. 2010. As e 
increases, the initial F value tends to increase (giving the characteristic exponential 
decay of the purity rate due to decoherence) while the amplitude of the oscillations 
seem to decrease approaching the value of the classical Lyapunov exponent. In fig- 
ured] (b) the decay rate of the purity Fj, which corresponds to the BE for E = 0). 
The coloured points correspond to the curves - for different e values - in figure [1] 
(a) . For the GDM we observe saturation of F^ at A as is expected. In the inset we 
assess the sum law F ^ F^ -I- Fs for the BE. There we plot F — F^ as a function 
of E/ft: the expected behaviour - all curves collapsing into the one corresponding 
to Fe - is only observed for values of e < 0.0035 corresponding to F^ < 0.5. For 
e = 0.0035 (0 symbols) we see that the sum law breaks up around E//i « 0.75. For 
e > 0.0035 the sum law is no longer valid. 



The next environment model that we considered is the generalised depolarising 
channel (DC). Although - as we shall see - in phase space it is somehow an ex- 
tremely non-local noise, its importance lies in that it is one of simplest and best 
known noise channels in quantum information formalism (Nielsen & Chuang 2000). 
The action of the DC for one qubit (A'' = 2) is simple: with probability (1 — e) it 
does nothing, and with probability e it 'depolarises' it, meaning that it leaves it in a 
completely mixed state. This is done by applying every possible Pauli matrix on the 
state. For an N dimensional system, and a torus phase space it can be generalised 
as follows (Aolita et al. 2004) 



that is, with probability (1 — e) it leaves the state unchanged, while with probability 
e it applies every possible translation in phase space, with equal weight e/N^. So 
contrary to the GDM where the incoherent sum over displaced states took place 
between states lying effectively close - due to the Gaussian weight-, for the DC the 
incoherent sum is over all states, close or apart. It is in this sense that we say this 
model is highly non-local. 

In figure [2] (a) we show the BE decay rate F as function of perturbation param- 
eter E for the perturbed cat map a = b ~ 2, N = 800 and k = 0.01 in the presence 
of DC noise model for distinct values of e. Again, here the red line with x symbols 
is Fs of the LE. For smaller e the curves look like essentially the same curve shifted 
upwards. There is no evident saturation at the Lyapunov exponent. For a larger e 



(&) Generalised depolarising channel 




(3.8) 
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Figure 2. (a) Decay rate F of the LB as a function of the rescaled strength of the 
perturbation E/ft for a DC environment. The map is the quantum version of the perturbed 
cat [equation (13. ip ] with a — b = 2. Averages were done over — 10 initial states. Other 
parameters are: k = 0.001, iV = 800. (X) e = (LE), { O) e = 0.1, (A) e = 0.22, (O) 
e = 0.40, (V) e — 0.7. The horizontal dashed (in (a) and (b)) lines correspond to the 
Lyapunov exponents of the corresponding maps A = ln[3 + 2-\/2] ~ 1.76275; (b) The decay 
rate Fe of the purity as a function of the perturbation parameter e. The points correspond 
to the initial values of the curves in (a) (inset) Decay rate F — Fe as a function of the 
rescaled strength of the perturbation S/ft. 



the BE oscillations tend to disappear and the growth is somehow linear with no 
apparent saturation. On figure [2] (b) we show the decay rate of the purity as a 
function of e and the coloured points mark the initial values of the curves on the 
top. Initially Fg grows linearly. As it is expected (Casabone et al. 2010), there is 
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Figure 3. (a) Decay rate F of the LB as a function of the rescaled strength of the per- 
turbation E/fi for a LDM environment. The map is the quantum version of the perturbed 
cat [equation (|3.ip ] with a = 6 = 2. Averages were done over Us — 10 initial states Other 
parameters are: k = 0.001, iV = 800. (X) e = (LE), ( □) e = 0.001, (A) e = 0.002, (O) 
e = 0.005, (V) e = 0.01. The horizontal dashed (in (a) and (b)) hues correspond to the 
Lyapunov exponents of the corresponding map A — ln[3 + 2\/2] ~ 1.76275; (b) The decay 
rate Fj of the purity as a function of the perturbation parameter e. The points correspond 
to the initial values of the curves in (a), (inset) Decay rate F — F^ as a function of the 
rescaled strength of the perturbation E/ft. 



no parameter independent regime for the DC observed, neither for F nor for F^. In 
the inset of figure [5] we show the decay rate of F — F^. We can see the lines collapse 
to the curve corresponding to Fs (red with x symbols) for the LE for a sizeable 
interval of E/fi, and up to values of F^ w 1. From the work by Casabone et al. (2010) 
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we know that the decay rate of purity as a function of e is = 2e, for small e. It 
is simple to show that in the interval of epsilon where this is valid holds the sum 
law Fs w r — Fe also holds. Here this is true up to values e < 0.4 (see also figure 2 
in Casabone et al. 2010) correspoding to F^ < 1. 



(c) Lorentzian decoherence 

Finally we consider a model which is more local than the DC but which unlike 
the GDM has polynomially decaying tails for c^{q,p). The motivation for using this 
model arose in the work by Casabone et al. 2010 when comparing the universalities 
of the purity and the LE. We take Ce{q,p) a sum of Lorentzians 



-1 ^ eiv 

Ce{q,p)^—r y 7 5 ^ (3-9) 



j,k= 

with A the proper normahsation for J2q p Ce(<Z,J') = 1 and q,p = 0, . . . ,N — 1. The 
sum is done to account for the periodicity of the torus (theoretically a; — > oo, prac- 
tically X is an integer much larger than 1). We call this model Lorentz decoherence 
model (LDM). Equation (|2.5p with Ce{q,p) given by equation (|3.9p defines a ran- 
dom process with Lorentzian weight that can be related to superdiffusion by Levy 
flights. The effect of heavy tails in decoherence is also explored in e.g the work of 
Schomerus & Lutz (2007). 

In figure we show the decay rate F of the BE as function of perturbation 
parameter S for the perturbed cat map a = b — 2, N — 800 and k — 0.01 in the 
presence of LDM for different e values. Again we see that for small e the curves 
look like a shift of one another - although less so than for the DC model- and then 
for large values of e the oscillations are destroyed and the growth of F is linear, like 
for the DC. On figure EJb) the decay rate F^ of the purity with the initial points 
of the curves on the top superimposed. The initial growth of Fj is quadratic with 
e as was shown in Casabone et al. (2010). It can also be clearly observed that in 
neither figure there is a parameter independent -Lyapunov - regime. In the inset 
of figure [3] we show the decay rate F — F^. The sum law Fs ~ F — Fj holds for an 
interval of S/fi, of up to w 1.5 (similar to the DC case) but it seems to break 
up a little bit earlier in the values of F^. Notice that in the inset of figure [Sjb), the 
line corresponding to the circles (F^ « 1) separates from the others at E/Zi « 0.75. 



4. Conclusions 

Summarising we have studied the BE for quantum chaotic maps with three dif- 
ferent types of decoherence. The BE complements the original idea of the LE in 
that it considers the presence of an environment yielding it appropriate for the un- 
derstanding realistic experiments. We have done extensive numerical calculations 
for a wide range of values of the perturbation of the map and the strength of the 
decoherence superoperator and we have focused on the decay rate of the BE in 
the regime where it decays exponentially. Other than providing a 'visual landscape' 
of the decay rate F of the BE our calculations enable a qualitative and quantita- 
tive analysis of the universal regimes found in the literature. We found that the 
more realistic diffusion model (GDM) correctly retrieves the Lyapunov behaviour 
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for large enough values of e. However, for this same case the sum law F « Fs + 
breaks up for relatively small values of e. We infer that this problem is related 
to the geometry of phase space (similar non-universal behaviour is found for the 
purity in the work of Casabone et al. 2010). On the contrary, the two other cases 
considered satisfy the sum law rather well. These two models have in common the 
slow decaying tails of the kernel c^{q,p), which means that the dccohercncc model 
acts non-locally in phase space. Furthermore these two models fail to exhibit the 
parameter independent Lyapunov regime. 

We have used quantum maps as generic chaotic systems and three very different 
decoherence models. We can thus conclude that non-generic behaviour is to be 
expected in echo experiments with arbitrary types environment. 

The authors acknowledge financial support from CONICET (PIP-6137), UBACyT (X237) 
and ANPCyT. D.A.W. and I. G.-M. are researchers of CONICET. 
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